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Abstract 

We experimentally demonstrate a sigmoidal variation of the composition profile across semi- 
conductor heterointerfaces. The wide range of material systems (III- arsenides, III-antimonides, 
III-V quaternary compounds, Ill-nitrides) exhibiting such a profile suggests a universal behavior. 
We show that sigmoidal profiles emerge from a simple model of cooperative growth mediated by 
two-dimensional island formation, wherein cooperative effects are described by a specific functional 
dependence of the sticking coefficient on the surface coverage. Experimental results confirm that, 
except in the very early stages, island growth prevails over nucleation as the mechanism governing 
the interface development and ultimately determines the sigmoidal shape of the chemical profile in 
these two-dimensional grown layers. In agreement with our experimental findings, the model also 
predicts a minimum value of the interfacial width, with the minimum attainable value depending 
on the chemical identity of the species. 

PACS numbers: 68.35.Fx, 68.43.Fg, 68.55.A-, 81.05.Ea, 81.15.Hi 
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A central goal of modern materials physics is the control of interfaces down to the atomic 
level. In particular, the behavior of layered materials depends on the atomic-scale structural 
roughness and chemical mixing across the interface Although abrupt interfaces between 
conventional semiconductors (such as III-V compounds) are fabricated and element profiles 
across these interfaces are obtained with atomic resolution, the relation between the layer 
growth processes and the parameters governing the interface formation and evolution is 
not satisfactorily understood. In this respect, there is an ongoing discussion about how 
interfaces can be quantitatively described on the basis of a growth model and whether there 
is a minimum interface width. 

Recently Hulko et al. {2], [J and Luna et al. (4 -f| have shown empirically that experi- 
mental concentration profiles in III-V two-dimensional (2D) heterostructures, e.g. quantum 
wells (QW) grown by molecular beam epitaxy (MBE), can be accurately reproduced by a 
sigmoidal function of the form x(z) = xq/[1 + exp(— z/L)]. Here, xq denotes the nominal 
mole fraction of one of the species, z is the position across the interface along the growth 
direction, and L is the parameter quantifying the interface width (L is proportional to the 
widely reported length W, over which the concentration changes from 10% to 90% of its 
plateau value). Moreover, the accuracy of the sigmoidal fitting seems to 
the experimental technique used to obtain the element distribution 
terestingly, of the compound semiconductor. In this letter, we show that a sigmoidal profile 
emerges from a simple model of cooperative growth with 2D island formation. Furthermore, 
the use of a generalized sigmoidal expression gives a reliable and systematic quantification of 
the chemical interface. It sheds light on basic aspects of the early stages of heteroepitaxial 
growth, and permits to find a correlation between the profile and the interface structure in 
morphologically perfect 2D epitaxial layers which have been grown in the 2D Frank- van-der- 
Merwe mode (FM), not necessarily by MBE. 

Experiments show that L depends strongly on the combination of materials on both sides 
of the interface, but that similar heterostructures feature the same value of L independently 
of their optimized growth methods [8( or substrate temperatures T s 9]. It could be argued 
that different T s might modify diffusion and intermixing processes at the interface and 
therefore account for the different values of L. This argument, however, does not explain 
the different values of L for dissimilar heterostructures despite their same T s or, conversely, 
that the same heterostructure features identical L when grown at different T s . Therefore, 
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the material interface property relation seems to be more complex than a mere increase in 
intermixing with increasing T s . Furthermore, the results suggest that there is a material- 
specific limitation in the interface width. The wide range of materials exhibiting a sigmoidal 
profile suggests a universal behavior determined by fundamental processes occurring during 
growth. 

A sigmoidal growth has been reported to occur in biological, geological and chemical 
processes with cooperative effects jlOj, wherein the binding of one atom or molecule affects 
the binding of the subsequent atoms or molecules. The experimental observation of sig- 
moidal profiles at semiconductor heterointerfaces suggests that similar phenomena occur in 
materials science. We infer that the species (atoms or molecules, hereafter we will use both 
terms indistinctly) which are involved in the growth and interface formation process form a 
strong cooperative system, where the binding affinity (quantified by the sticking coefficient 
s) changes with the amount of atoms that have already been bound. More concretely, we 
assume that the rate of adsorption of A atoms (let us consider the growth of a ternary 
semiconductor alloy A X0 B\- X0 C on top of a binary alloy BC) is proportional both to the 
surface concentration of chemisorbed A atoms a and to the surface concentration of free A 
sites (er — where o~ is the total surface concentration of A sites: 

^= S -lla(t)(a -a(t)). (1) 
at o"o 

Here si is a parameter that depends on the material system (we will see later that it cor- 
responds to the initial island growth sticking coefficient) and J is the impingement rate on 
an area of the substrate corresponding to one adsorption site 11, [l2]. If we assume that, 



after fast local rearrangements not described by a simple kinetic equation, the thickness z 
is proportional to the time variable t, then the variation of the surface coverage of A atoms 
6 = a /(To is described by the differential equation: 

§ = ij£, W( l -,(,,), (2) 

where the constant r is the mean growth rate of the ternary alloy. The solutions of Eq. (J2J) 
with initial conditions between zero and one are sigmoidal functions. Equations (Q]) and (J2J) 
feature a surface-coverage dependent sticking coefficient s(9) = Si#(l — 6). Incidentally, we 
mention that the experimentally measured value of the mole fraction x is simply proportional 

to e. 
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(a) Nucleation 



(b) Island mediated growth 
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FIG. 1. (Color online) Schematic location of the chemisorbed atoms on a surface when (a) nu- 
cleation and (b) island growth dominates. The island-mediated growth proceeds via Eden cluster 
formation: in each growth step, empty perimeter sites next-neighbored to a seed (or to an occupied 
site) become occupied. 



Comprehensive investigations on the MBE growth of the II- VI compound CdTe showed 
that the sticking coefficients for both the Cd and the Te atoms were not constant jl3 |. 
Based on a remarkably good fitting of their data, Litz et al. suggested the existence of a 
precursor state during MBE growth, in which atoms are physisorbed to the growing sur- 
face in the manner of Kisliuk's theory of precursor-mediated adsorption . In his theory, 
Kisliuk distinguished two types of physisorbed molecules, intrinsic precursor molecules (ipm) 
and extrinsic precursor molecules (epm), which are physisorbed above vacant and occupied 
chemisorption sites, respectively, albeit he assumed that adsorbed molecules are uncorre- 
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14|. Ensuing computer simulations revealed that molecules adsorbed via an epm 
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mechanism tend to aggregate into 2D clusters or islands [15] . Becker et al. 
Kisliuk's model to include these lateral positional correlations. In their model, Becker et al. 
addressed chemisorption as a combination of two mechanisms, each corresponding to the 
two types of precursor states introduced by Kisliuk: the ipm are responsible for a seeding 
process in which a molecule is directly chemisorbed above a vacant site and serves as a 
nucleus for the formation of an island, while epm are responsible for island growth by initial 
physisorption above an occupied chemisorption site followed by lateral jumps, until a vacant 
chemisorption site is reached. In this latter process, attractive adsorbate-adsorbate inter- 
actions enhance the precursor binding potential near the island edges, so that the epm will 
most likely become chemisorbed next to an already chemisorbed atom at the island edge, 
generating 2D compact Eden-type clusters 
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17j . as it is illustrated in Fig. [TJ Further- 



more, Becker et al. proposed the additive contributions of the (ipm) nucleation and (epm) 



island growth terms to the sticking coefficient: for the nucleation process, they assumed a 
Langmuir form where the sticking coefficient is proportional to the fraction of empty sites, 
i.e., to Sq(1— 9) where Sq is the initial sticking coefficient [18(, while for the island growth term 



they assumed dependencies such as yl — 9, which do not lead to pure sigmoidal profiles. 

To gain insight into the phenomena underlying the interface formation, we have gener- 
alized Eq. (j2J) by adding the Langmuir dependence of the seeding process (which describes 
independent nucleation sites that adsorb no more than one adsorbate and leads to a random 
spatial distribution of nucleation centers on the surface) to our growth term (jSJ). The latter 
describes the attractive adsorbate-adsorbate interaction previously mentioned through the 
9(1 — 9) dependence. Thus, Eq. (J2J) is extended to: 

HQ T 

- = -[so(l-9(z)) + Sl 9(z)(l-9(z))}, (3) 
dz r 

or, equivalently, 

df) 1 

Equations (j3J) and (j4]) feature the generalized sticking coefficient 

s(9) = s (l-9) + Sl 9(l-9), (5) 

the fraction 9q = So/ s i denotes the ratio of the initial sticking coefficients from the nucleation 
and island growth terms, respectively, and 

l = _lH_^ (6) 

SO + S! 

We remark that this equation gives the aimed connection between the interfacial profile, the 
growth process, and the material system. In particular, it gives the quantitative relation be- 
tween the interface width L, the kinetic parameters r and J, and the microscopic parameters 
So and S\. Furthermore, it refines the definition of an abrupt interface. Since So + s% < 2, a 
non-vanishing interface width is predicted for any semiconductor heterointerface. 
The general solution of Eq. (j4j) can be written in the form 

e(z) = i - - l ] 6 \ IT ) (7) 

V ' I-)- e {z-z )/L ' V > 

where z Q is an integration constant that merely shifts the z coordinate. Information on the 
main processes governing the growth and interface formation can be gained from inspection 
of Eq. ([7]). In Fig. |2] we plot the profiles obtained from Eq. (JTJ) with 9(0) = [which 
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FIG. 2. Variation of the profiles given by Eq. ([7|) as 9q = sq/si ranges from 10 -2 to 10 2 . The limit 
#o 3> 1 corresponds to nucleation, while the limit 6q <C 1 corresponds to island growth. 

corresponds to zq = — Lm(#o)] and 6q ranging from 1CT 2 to 10 2 . The asymmetry of the 
profiles is a quantitative measure of the relative weight of the nucleation and the island 
growth terms: in the limit 6 3> 1, i.e., So 3> Si, where nucleation is the only mechanism, 
the profiles are very asymmetric and tend to 9{z) = 1 — exp(— z/L); the asymmetry of the 
profiles decreases with 8 , and at 6 = 1, i.e., when s = si, an inflection point appears 
in the profile; finally, in the limit 6q <C 1 where island growth dominates, we recover the 
pure, symmetric sigmoidal profile 9(z) = 1/[1 + exp(—z/L)]. Note, however, that if so = 
the solution of the differential equation (j2j) with 6(0) = is identically zero, which shows 
that nucleation must always exist. However, in this case, 2D island growth occurs from 
the very early stages of the interface formation, since very few but already existing nuclei 
are sufficient to promote further atom attachment at their rims (Eden growth, cf. Fig. [1]). 
Obviously this process is initially less efficient than seeding, and initially it takes more time 
to increase the composition profile to significative values. 

We have analyzed experimental element profiles using the functional form in Eq. (J7|) and 
derived 6q an d L for different III-V semiconductor heterostructures. We have investigated 
more than 60 samples comprising 2D morphologically perfect Ill-arsenides, III-antimonides, 
III-V metastable quaternary compounds and Ill-nitrides heterostructures, which were fab- 
ricated using different methods. High-resolution and analytical transmission electron mi- 
croscopy (TEM) techniques have been used to quantitatively determine the change in sto- 
ichiometry across the interfaces. Representative results are summarized in Table [B As ob- 
served, in all cases, regardless of the materials system and of the growth method, the largest 
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TABLE I. Parameters obtained from nonlinear fittings of the experimental profiles of 2D morpho- 
logically perfect III-V semiconductor heterostructures to the theoretical profile of Eq. (JT]). Growth 
methods: molecular beam epitaxy (MBE), solid phase epitaxy (SPE), metal-organic chemical va- 
por deposition (MOCVD). References are given for data taken from the bibliography. As a measure 



of the overall quality of the fit we include in the last column the coefficient of determination R 2 . 



Material 


Ref. 


Growth method 


x (%) 


8o 




L (ML) 


R 2 


Al xg Gai - X[t As 1 GaAs 




MBE 


29.46 ± 0.02 


(1.7 ±0.8) x 10" 


3 


1.65 ±0.02 


0.9997 


In^oGai-xoAs/GaAs 




MBE 


39.8 ±0.06 


(5.5 ± 0.7) x 10~ 


3 


1.31 ±0.01 


0.9996 


In^ Gai-^As/ GaAs 




SPE 


34.2 ±0.02 


(2.3 ± 1.8) x 10" 


3 


1.61 ±0.03 


0.9991 


In^oGai-^Q As/GaSb 


[19] 


MBE 


49.2 ±2 






1.0 ±0.1 




Ini-^GaajoSb/InAs 


[20] 


MBE 


81 ±3 


(1 ± 1) x 10~ 2 




1.3 ±0.1 


0.9970 


(In, Ga)(N, As)/GaAs 




MBE 


(In) 32.6 ±0.2 


(2.1 ± 2.2) x 10" 


3 


1.62 ±0.03 


0.9990 








(N) 5.7 ±0.1 


(1.5 ± 0.7) x 10~ 


2 


1.3 ±0.1 


0.9892 


(In, Ga)(N, As)/GaAs 




MBE 


(In) 34.4 ±0.1 


(2.9 ± 2.5) x 10" 


3 


1.45 ±0.03 


0.9996 








(N) 4.33 ± 0.06 


(6.7 ±5.4) x 10" 


3 


1.0 ±0.1 


0.9988 


In^Gai.^N/GaN 




MOCVD 


23.9 ±0.2 


(1.8 ± 0.4) x 10- 


2 


1.5 ±0.1 


0.9986 


In^Gai.^N/GaN 


[22] 


MOCVD 


7.15 ±0.02 


(1 ± 1) x 10- 3 




(0.177 ±0.002) nm 


0.9995 



values for 9q are 9q ~ 2 x 10~ 2 , which after due account of the accuracy, are consistent with 
a pure sigmoidal appearance of the experimental curves. Additionally, L varies with the 
combination of materials at both sides of the interface. As an example of these data and of 
the resulting fits, we show in Fig. E)[a) an experimental profile across an (Al, Ga)As/GaAs 
interface with L = 1.6 ML, while Figs. [3](b) and [3](c) correspond to the In and N dis- 
tributions, respectively, across an (In, Ga)(As, N)/GaAs QW. In this case, the fit yields 
L(In) = 1.6 ML and L(N) = 1.3 ML for the (In, Ga) (As, N) /GaAs interface, while for the 
GaAs/(In,Ga)(As,N) interface we obtain L(In) = 1.7 ML and L{N) = 1.3 ML (cf. Table CD). 
In all our experiments, we have not found neither significant asymmetric distribution with 
an exponential-like behavior nor curves reflecting ^ « 1. Furthermore, the experimental 
data are remarkably well reproduced by sigmoidal profiles, like those represented in Fig. |3j 
We thus conclude that, in spite of the ubiquitous nucleation term, the epm-mediated 2D 
island growth mechanism is the primary process at the forming interfaces and determines 
the general shape of the transition zone at the heterointerface. Note that, although the exact 



adsorption processes occurring in the epitaxy of semiconductors are generally unknown 23) , 



it is very well established by reflection high-energy electron diffraction or scanning tunneling 
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FIG. 3. (Color online) Experimental composition profile (open symbols) and sigmoidal function fits 
(solid line) across (a) an interface of (Al,Ga)As/GaAs. (b) and (c) represent the In and N profiles, 
respectively, across the interfaces in the QW structure (In,Ga)(As,N)/GaAs. The experimental 
relative error of the data points is typically 1.5% for Al and In and 4.6% for N. (The growth 
direction goes from left to right). 



microscopy measurements that in 
in an island-mediated 2D process 



;he FM growth mode the completion of the layers occurs 
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Furthermore, our comparative analysis of several III-V heterostructures reveals a limita- 
tion in the interface width, which is specific to each material system, in agreement with the 
predicted limit in the value of L. As shown above, the interface width L contains information 
both on the growth kinetics and on the adsorbate-surface relation (through sq and s\). In 
general, the magnitude of the sticking coefficient and the adsorption dynamics are ultimately 
determined by the exact details of the potential energy surface (PES) that describes the par- 
ticular molecular-surface interaction and is specific to a given molecule-surface system 18]. 
Our experimental findings suggest that, for samples grown under optimized conditions [8j, 
the minimum attainable interface width is not kinetically driven, but thermodynamically 
driven through the adsorbate-surface interaction. The experimental results support this 
hypothesis. The width L at the A a . B 1 _ x . () C/BC interface would be determined ultimately 
by the A a d SO rbate - BCsurfaco PES, which depends on the chemical identity of the surface 
and adsorbing atoms. Hence, the different values of L for (Al, Ga)As/GaAs compared to 
(In, Ga)As/ GaAs (cf. Table I) may arise from the different PES in the Al a d SO rbate - GaAs sur faco 
and In adsorbate - GaAs surface systems, respectively; and, similarly, the In adsorbate - GaAs surfaco 
interaction will differ from the In adsorbate - GaSb surface one, therefore leading to distinct val- 
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ues of L in the (In, Ga)As/GaAs and (In, Ga)As/GaSb material systems, as experimentally 
observed. Finally, the intriguing identical L for high-quality (In, Ga) (As, N) / GaAs QWs re- 
gardless of the growth procedure or laboratory [9J] is also explained because the PES would 
be similar in all cases. Equation (jH]) also predicts a dependence on the substrate orientation 
because the PES (and hence Sq and S\) depends on the crystallographic orientation of the 
surface above which growth proceeds. 

Finally, we would like to point out that there is a close relation between Eq. ([3]) and a 
similar equation obtained in the context of silicon oxidation 25], where a sigmoidal behav- 
ior was also reported. Although seemingly dealing with a different problem, it has been 
recently shown that the initial stages of Si(OOl) oxidation are determined by a combina- 
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26| in close analogy with our 



tion of a Langmuir-type adsorption and 2D island growth 
study. Moreover, sigmoidal laws have been reported in the context of nanocrystals growth 
in solution [27], where recent investigations reveal a complex growth dynamics involving a 



2D nucleation and growth process via atomic attachment from the solution phase [28|, 129 ]. 

We conclude that a sigmoidal variation of the composition across an interface (or, in 
general, of any order parameter defining a system) directly arises from this specific growth 
dynamics characterized by strong cooperative interactions with 2D island formation and, 
therefore, it has a general character and is not restricted to MBE or MOCVD grown semi- 
conductor heterointerfaces. We show that the prevalence of precursor-mediated 2D island 
growth over nucleation is the primary mechanism governing the interface development in 
2D semiconductor heterostructures and responsible for the sigmoidal shape of the chemical 
profile at the interface. Moreover, while kinetics factors confirm as valuable control param- 
eters for adjusting the interface width during the fabrication process, there is a minimum 
attainable interfacial width, which is dictated by the molecule-surface interaction potential 
and is exclusive for each material system considered, since the latter depends on the chemical 
identity of the species involved. 
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